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Saturated magnetic field amplification at supernova shocks 
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ABSTRACT 

Cosmic-ray streaming instabilities at supernova shocks are discussed in the quasilinear 
diffusion formalism which takes into account the feedback effect of wave growth on the 
cosmic ray streaming motion. In particular, the nonresonant instability that leads to 
magnetic field amplification in the short wavelength regime is considered. The linear 
growth rate is calculated using kinetic theory for a streaming distribution. We show 
that the nonresonant instability is actually driven by a compensating current in the 
background plasma. The nonresonant instability can develop into a nonlinear regime 
generating turbulence. The saturation of the amplified magnetic fields due to particle 
diffusion in the turbulence is derived analytically. It is shown that the evolution of 
parallel and perpendicular cosmic-ray pressures is predominantly determined by non- 
resonant diffusion. However, the saturation is determined by resonant diffusion which 
tends to reduce the streaming motion through pitch angle scattering. The saturated 
level can exceed the mean background magnetic field. 
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1 INTRODUCTION 

Diffusive shock acceleration (DSA) is regarded as the pre- 
ferred mech anism for the acceleration of the Galactic cosmic 
rays (CRs) (|Drurvlll983l ). It is commonly believed that SN 
shocks are responsible for acceleration of high energy CRs 
at least up to the 'knee' (~ 4 x 10 eV) of the CR spectrum 
(jHillas! [20061 ) . In the DSA model, particles gain energy by 
bouncing back and forth across the shock. Although parti- 
cles only gain a small amount of energy in each crossing, 
they can be accelerated to very high energy through many 
crossings provided they can be trapped in the acceleration 
region for a sufficiently long time. The acceleration is rather 
efficient and naturally leads to a power-law energy distri- 
bution. Despite these advantages, there have been two long 
standing problems. First, the standard DSA theory pr edicts 
the maximum energy well below the 'knee' (Bell 2004). The 
maximum energy achievable is limited by both the Bohm 
approximation, in which the particle's gyroradius must not 
exceed the mean free path to scattering, and the condition 
that the gyroradius is smaller than the shock's width (other- 
wise the particle would escape from the acceler ation region) 
|Bellll2004 IZirakashvili. Ptuskin fc Volkll200sl ). For a mag- 
netic field of order 10 -10 T and typical parameters of SN 
shocks, one has the maximum energy ~ 5 x 10 14 eV. Second, 
turbulence is required for effective scattering in the acceler- 
ation region so that the particles can be trapped (|Skillind 
1 19751 ; lBel3ll978l : lLagage fc Cesarskvlll983al lbl). So far, in the 
standard DSA theory, one generally postulates that Alfven 



turbulence can be generated by CRs themselves through res- 
onant interactions, even though in practice the growth of 
Alfven waves due to CRs is known to be ineffective. 

There is growing interest in the possibility that the mag- 
netic field at the shock may be amplified due to CR induced 
instability in the nonre s onant regi me dMcKenzie fc Volkl 
1 19821 ; iBell fc Lucekll200ll ; lBellll2004l ). iBelll i|2004h showed a 
nonresonant form of the instability driven by a CR current 
can outgrow the familiar resonant form, leading to mag- 
netic amplification. A strong magnetic field at shocks re- 
duces the gyroradius and this raises the maximum energy 
to which particles can be accelerated. X-ray observations 
of young SNRs suggest that the magnetic field strength 
near the shocks is much higher than that in the ISM and 
is a s trong function of the shock speed (IVink fc Laming! 



spe 

120031 ; I Yolk. Ksenofontov fc Berezhkol [20ol ) . The presence 
of strong magnetic fields may explain the lack of strong 
TeV gamma-ray fluxes, as shown fro m HESS observations 
l|Volk. Ksenofontov fc Berezhkol I2OO8T ) . The nonthermal X- 
rays are well described by synchrotron spectra. A stronger 
magnetic field in the emission region implies that a lower 
electron energy is required. This would lead to a lower TeV 
gamma flux from inverse Compton scattering. 

The nonresonant growth of Alfven waves due to CR 
streaming has been discussed in both the MHD formalism 
llBell fc Lucekll200ll ; lBen1l2004l ; IZirakashvili. Ptuskin fc Volkl 
20081) and the k i netic f ormalism dMelrosd 120051 ; 



Reville. Kirk fc Duffvl 120061 ; lAmato fc Blasil 120081 ). In 



both formalisms, the instability is shown to exist in the 
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linear regime. To test the magnetic field amplification model 
against observations one needs to determine the saturated 
magnetic field accurately. The saturation of the instability 
cannot be determined in linear theory as the reaction of 
the instability on the CR momemtum distribution is not 
automatically included in the linear calculation. Although 
there are numerical simulations of the instability that 
extend to the nonlin ear regime, different saturation levels 
have b een predicted (|Niemiec. Pohl. Stroman. fc Nishikawal 



120081 ; iRiauelme fc Spitkovskvl 120081 ). In this paper, the 
instability is discussed in the quasilinear formalism in which 
the reaction of the instability on the CR distribution can 
be included self-consistently. So, in this formalism one can 
estimate the saturation analytically, with both nonresonant 
and resonant diffusion processes considered. The treatment 
of the nonresonant diffusion present ed here is simil ar to 
that used for the firehose instability (|Davidsonlll972l ). We 
emphasize the major difference between the CR streaming 
instability and the firehose instability. The former is caused 
by streaming motion and the latter is due to a pressure 
anisotropy with excess of parallel pressure over the perpen- 
dicular pressure (with respect to the mean magnetic field). 
To some extent, the streaming instability also resembles the 
Weibel instability — a nonresonant, purely growin g mode 
drive n by anisotropy in the particle distribution l|Weibell 
1 19591 ). 

In Sec 2 the kinetic theory of CR streaming instabilities 
is discussed with emphasis on the nonresonant instability. 
Quasilinear diffusion driven by the nonresonant instability 
is discussed in both the short and long wavelength approxi- 
mations in the nonresonant regime in Sec 3 and in the reso- 
nant regime in Sec 4. Application to SN shocks is discussed 
in Sec 5. 
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For the standard DSA one has p = 2 (|Belfll978l : rDrurvlll983r i 
and when the nonlinear effe ct on DSA is i ncluded, p deviates 
from this canonical value l|Eichlerl 1 1984 ). The distrib ution 
with a = 1 corresponds to that used in Melrose! (|l986h . The 
distribution |T]) implies 



n-CR = 4-7T 



= 2tt 



U±dUA 



sin ada 



u du f(u cos a, u sin a) 



(3) 



where cos a = uu/u. In the second expression in (|3}, one 
chooses (u, a) in place of (uii,u_l) as independent variables. 
It can be verified that averaging the parallel velocity v\\ 
over the distribution (JXJ gives the streaming velocity vcr, 
which is independent of the choice of the parameter a. The 
CR current is then given by Jen = qncRV CR . The pres- 
ence of streaming CRs affects the background plasma in two 
ways, due to their charge density and their current density, 
respectively. The background plasma must have a charge 
density and a current density that are equal and opposite 
to those of the CRs. This requires that the electrons and 
ions (assumed to be protons) have different charge densi- 
ties, n e 7^ rip, and that they move relative to each other 
with streaming velocities v e 7^ v p , which are assumed to be 
along the gui ding magnetic fiel d. The neutralization condi- 
tions require |Achterberdll983l ) 



e(n e 



qn c 



e(n e v e ~n p v p ) = qn CB v CR . 



(4) 



These properties of the background plasma drive the non- 
resonant instability attributed to the CRs. 



2 COSMIC-RAY STREAMING INSTABILITIES 

We outline the kinetic theory of CR streaming instabili- 
ties including both the usual resonant instability and the 
nonresonant instability and focus particularly on the lat- 
ter. Our treatment builds on other recent discussions of lin- 
ear kinetic theory of the CR-induced nonresonant instability 
l|Reville. Kirk fc Duffvll2006l ; lAmato fc Blasill2008l ). For con- 
venience we assume a single species of CRs with charge q and 
mass m. 



2.1 CR streaming motion 

To model CR streaming at velocity vcr along the mean 
magnetic field we consider a class of streaming distributions 
in momentum space given by 



/(u||,«_l) = n C R 



(a , + i)!!S»(IilY 

V \ u J 
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where Uu and u± are the nondimensional momenta (normal- 
ized by mc) parallel and perpendicular to the mean magnetic 
field, respectively, a is an integer 1, ucr is the CR num- 
ber density, v is the CR's velocity written as a function of 
My and u±, u = (u\ + ity) 1 / 2 , and 



2.2 Dispersion relation 

A formal procedure to derive the dispersion relation involves 
separating the plasma response tensor into that for a back- 
ground plasma, denoted by Kij, plus that for the CR com- 
ponent, denoted by Aify. The background plasma can be 
regarded as a cold, magnetized plasma, while the CR com- 
ponent is described by the distribution ([T]). Assume that 
the gyrofrequency of CRs is SI — \q\B/m, where B is the 
mean background magnetic field. A useful approximation is 
k±v±/£l <C 1, where k± is the perpendicular wave number, 
tl = fi/7 and 7 = (1 + u 2 ) 1 / 2 is the Lorentz factor of CRs. 
The approximation implies that in the response tensor, only 
the first gyroharmonics terms are important. For the back- 
ground plasma, one assumes v\ <C c 2 and the low-frequency 
approximations, u) tti, ui <C |fc||f||| and to <g Q, where Qi 
is the gyrofrequency of ions in the background plasma. The 
background magnetic field is assumed to be along the 3-axis. 
Since K33 oc Sli/co 2 can be set to 00, only the 2x2 compo- 
nents of the response tensor are relevant. For the background 
plasma, these components can be written as 



Ku — K22 = —5-, 



-K21 



2n t 



i p fc||(n e u e 
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turii 
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where nil = knc/u), va = fli/uipi is the Alfven speed, cj p i is 
the ion plasma frequency of the background plasma. Using 
the neutralization conditions (UJ), © can be rewritten in the 
form 

HqJcr 



Kl2 = -K21 = 



• 2 (0) 

- l Vn\\X h , 



(0) 
Xh = 



2k«B 



(7) 



The CR components are 
AK U = AK22 = 



AK 12 = -AK 21 = -in\ X hM 



where 



Xa = -A I dxx 



■ 1 

Xh — A I dxx 

CT-1 



p+ 2ct - 3 (1 - X 2 )H(1 - x 2 ) 
1 + x 



(9) 



p+2a-3 



(1-x )ln 



l-x 



+ 2x 



4s 2 - 1 



A = 



iioJcR (4<7 2 -l)(p-l) 



4k«B 



6(l/wi, I/U2), b 



(10) 
(11) 

fi/|fci||c and H(x) the 



with {x-l,X2 
step function. 

Eq ([9]l and ()10|) correspond to the dissipative and reac- 
tive part respectively of the response tensor due to CRs. In 
the strong magnetic field limit x > 1, the dissipative part is 
zero. When x < 1, both dissipative and reactive terms con- 
tribute to the dispersion. The dissipative part is due to CRs 
in the resonance x = b/u ~ fi as discussed extensively in 
the context of the resonant instability (see discussion below) 
llKulsrud fc Pearcdll969| : iMeh-ose fc Wentzellll97d: ISkilli"n3 
Il975l ; iMcKenzie fc Volklll982l ; lLagage fc CesarsQ Il983al ) 
The dispersion relation takes the following form 



¥ + 1 + 2i Xa 



1/2 



, 2 2 



In the Xa — > and X n —* limit (xi must also be zero 



(12) 



(|12[) with the lower sign reproduces the usual Alfven mode 
dispersion u — \ku \va and with the upper sign the fast mode 
dispersion lo = kvA- 

2.3 Streaming instabilities 

Instability occurs when u> has an imaginary part of the ap- 
propriate sign. An imaginary part can arise in two different 
ways, which we refer to as resonant and nonresonant insta- 
bilities. Writing lo — > cu+iT in (|12[) and assuming r<w, the 
resonant instability is described by equating the small imag- 
inary terms, giving V — Xa k^v A /2uj. The value of Xa is de- 
termined by the CR distribution. A resonant instability can 
develop for both Alfven and fast-mode waves if streaming 
CRs are in resonance with the wave considered. For Alfven 
waves (lo — \k\\\v A ), the growth rate is 



(13) 



The polarization of both modes is generally elliptical and 
becomes approximately circular for nearly parallel propa- 



gation |fcx/fc||| <C 1. In each mode, waves of both sense of 
polarization can grow. 

On neglecting Xa, (|12[) with the upper sign becomes a 
real equation, because the quantity inside the square root is 
a sum of squares. For the negative sign (Alfven waves), lo 2 
can be negative, and one of the two solutions corresponds 
to an intrinsically growing wave. This is identified as the 
nonresonant instability. For |fc_L/fcii| 1, the condition for 
this nonresonant instability to occur is (xi°' +Xh) 2 > 1- The 
growth rate is 



V2 



(°) 1 ! 

X h +Xh\ 



- 1 



1/2 



fc 



\\\VA. 



(14) 



In principle, this can be satisfied either due to \xh\ > 1 or 
to \Xh°^\ > 1- F° r o = 1, the integral Xh given by (|10[1 can 
be calculated exactly, i.e., 



3^Jc«(6- l) , (/( (( 
4fc||Bx^ 



(l-z )ln 



1+x 



l-x 



+ 2x 



16fciiBxi 



/ 2 



Xl(l + Ij) - x 2 (l + x 2 ) 



1 + Xi 



1 — X\ 



1 I 2 
+ (X 2 



l) 2 ln 



l + x 2 



X2 



(15) 



The second equality in (|15p is obtained for p = 2. The con- 
dition \xh\ > 1 cannot be satisfied for x\ <C 1. Thus, only 
the second possibility ( X h ^ 1) ' s relevant. 
For 12 < 11 < 1, one finds 



Xh 



4k n B ' 



(16) 



so that the contribution from the CR current is much smaller 
than the contribution due to the compensating current in the 
background plasma. It follows that in treating the nonres- 
onant instability, one can neglect the direct contribution of 
the CRs, described by Xh, in comparison with the indirect 
contribution, described by xi°^ The growth rate (|14|) can 
be written in the following approximate form 



V2(xL 0) 



1/2 



t\\VA 



V2 



HqJcr 
kuB 



1/2 



Nonresonant instability requires 

\Jcr\ > 

Mo 



fc||«A-(17) 



(18) 



An explanation for why the direct contribution from the CR 
current is so much smaller than the indirect contribution 
from the background plasma is that in the limit feiir g ^> 1, 
most of CRs move rather rigidly. Their only role in this limit 
is to induce the compensating current in the background 
plasma that drives the instability. It is appropriate to point 
out here that a small fraction of CRs can satisfy the resonant 
condition (/1 = x) and that resonant instability due to these 
resonant CRs is insignificant compared to the nonresonant 
instability due to the compensating current. However, the 
resonant interactions between these CRs and waves with 
fcii 2> l/r g can lead to resonant diffusion that has feedback 
effects on CR streaming motion (cf. Sec. 4). 
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3 QUASILINEAR DIFFUSION 



1, 2, 3. The n ^ 2 terms can be avoided using 



To determine saturation we calculate the back reaction of 
wave growth on CRs in the quasilinear diff usion theory 
i|Shapiro fc Shevchenkol fl96l ; lDavidson|[l972l ). We assume 
that the nonresonant instability discussed in Sec 2 can de- 
velop well into the nonlinear regime, producing turbulence. 
The growing turbulence causes CRs to diffuse in momentum 
space and this in turn reduces the anisotropy, suppressing 
the instability. 

We calculate the particle diffusion in the weak turbu- 
lence approximation. The distribution of CRs in momentum 
space can be written as f(u,t) = F(u«,u±,t) + f^(u,t), 
where is a linear function of the electric and magnetic 
fields of the wave and F is the mean distribution, which is 
assumed to be axisymmetric with respect to the mean mag- 
netic field, B, and as the zeroth order approximation it can 
be taken to be ([T]). The response of particles to wave growth 
is determined by including quadratic terms of the electric 
and magnetic fields of waves. The diffusion equation for F 
can be obtained from the Vlasov equation (the derivation is 
outlined in Appendix). We assume k±r g <C 1 so that only 
the first gyroharmonics terms are relevant. This condition is 
generally satisfied for waves propagating nearly parallel to 
the mean background magnetic field. 



$2 = -* 1 + T7jP 

9$i _ it || <9$i 
chin b db ' 



d$ 2 _ 3u|j 9<£>i 
dun b db 



My d 1 ggi 
~bdbb~db 



(21) 
(22) 

(23) 



Multiplying (|19p by u| and u\ , respectively, and integrating 
them over 2-KU±du±duu divided by ncn, one obtains 



d(u\) 
dt 



d{u 2 ± ) 
dt 



where 



1 

27 



2 



(( 

ld_ld_ 
bdbbdb 



•4uf)$i) + 



2d_ 

bdb 



((2«i-«f)uf*i) 
(24) 
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27 



1 d 



2(^> + ^(«I*i(2«l-3«i)) 



13 13/4. 2 \ 



(25) 



3.1 Nonresonant diffusion equation 

For a nonresonant instability, it is relevant to consider the 
approximation that the resonant width is much larger than 
the growth rate, i.e., \kuvu — 3> F, where T > is the 
growth rate and for convenience the frequency is set to zero. 
The diffusion equation can be expanded in r/[ft||VM — £2/7|. 
Let 5Bk be the spatially Fourier-transformed magnetic fluc- 
tuations arising from the instability. The energy of mag- 
netic fluctuations in an elementary phase volume dfc/(2-7r) 3 
is \SB k \ 2 /2fj, . The expansion yields the following approxi- 
mate diffusion equation (Appendix A) 
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dt 
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6 2: 



ui+b 2 



b 2 y 



(19) 



(20) 



with b = B/B*, B, 
region considered. 



|fc|||mc/|g| and V the volume of the 



3.2 Perpendicular and parallel kinetic energy 

ft is of interest to examine the evolution of the parallel 
and perpendicular kinetic energy by evaluating the rate of 
change in the average u 2 and v?^ (over the mean distribution 
F). The averaging involves singular terms l/(w|| ± b) n with 



dk 2T \SB k \ 2 



(2tt): 



Bl 



(26) 



The ratio of the parallel to perpendicular pressures is de- 
termined by the same angular averages as the ratio of (u 2 ) 
to (wj_/2). It follows that the ratio of Eqs fl24} and 
effectively determines how the ratio of the parallel to per- 
pendicular pressures evolves. 



3.3 Long wavelength approximation 

In the long wavelength approximation 1/kn 3> r g , one has 
b ^> I liii I and 



$1 



(27) 



The rates of change in the average of ti» and u\ is given by 



d(u\) 

dt 
d{u\) 

dt 



_( 2 (u|>-(tti>)i, 



(28) 
(29) 



where r' has the same form as (|26p but with B, replaced by 
B. Thus, one has 1/r' = d(SB 2 / B 2 )/dt. One may compare 
the nonresonant diffusion to the firehose instability. Similar 
to the firehose instability, a growing wave causes the parallel 
kinetic energy to decrease and the perpendicular energy to 
increase with time. For the CR streaming distribution ([1]), 
one can show that the right-hand sides of (|28|) and (|29p are 
zero. In the long wavelength regime an instability can arise 
from both dissipative (resonant) and reactive (nonresonant) 
effects. As we focus on the nonresonant instability in the 
short wavelength regime, quasilinear diffusion in the long 
wavelength regime is not discussed further. 
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3.4 Short wavelength approximation 



The right-hand side of Eq (|24[) and (|25p can be evaluated in 
a general case in which the particle distribution is isotropic 
in pitch angles. The angular integration in (|24|l and (|25[) 
leads to 



d(uf) 
dt 

d{u\) 
dt 





1 - X 


- ( 2 + x In 






1+x 



-i/ | + 2a: 2 + x 3 In 



(30) 



(31) 



where one uses the following integral 



I 2n = u 2 / ^ 2n $id/i 



= 2 E 



1 -X 



1+X 



(32) 



The second expression in (|32[) is obtained by retaining the 
Cauchy principal value. Surprisingly, Eq (130 ^ and (|31|) have 
an opposite sign compared to their counterpart in the long 
wavelength regime (Eq 1281 and I29 [) . Although both the fire- 
hose instability and the instability discussed here are driven 
by anisotropy in the particle distribution in momentum 
space, there is an important difference between the two. The 
firehose instability is driven by anisotropy in kinetic energy 
with {u 2 } > (it 2 L )/2. By constrast, the streaming instabil- 
ity is due to the CR streaming motion ((uii) 7^ 0), which 
is also called streaming anisotropy, and the instability can 
develop even when the CR pressure is isotropic. With the 
specific choice of the distribution {l}, one can show that 
(it 2 ) = {u\)/2. Since the current in |T} is oc /i 2<T_1 , it does 
not contribute to the average (Y) if Y is proportional to an 
even power of fi. 

The instability can lead to pressure anisotropy with per- 
pendicular pressure increasing. In the approximation 1 < 1, 
one has d{u 2 )/dt ~ 2/r and d{u\) /dt ~ — 2/3r. This gives 



d 
di 



(2{ul)-{u 2 ± )) 



> 0. 



(33) 



One should emphasize here that development of such 
anisotropy in kinetic energy is driven by the streaming mo- 
tion of CRs. 

From (|19[). one may calculate the rate of change in the 
streaming speed due to nonresonant diffusion. One starts 
with a particle distribution with streaming motion similar to 
{TJ but without specifying the specific form of g(p)- Consider 
the case a = 1, i.e., the streaming component is oc /1. As in 
Sec 4.4, in the averaging, one can carry out the angular 
integration first. This gives 



-3^(1- ^1. + 
+/3 2 (6-7/3 2 /3 2 )/6 + 3/3 4 /, 
1 d 1 d 



-3/2 + (1 + 5/3 2 )/ 4 



+ 



x dx x dx 



P 2 h 



(1 



(34) 



where l2 n is given by (|32 
limit i<l. one has 



Since dv CK /dt = d{v^)/dt, in the 



dv r . 



3«, 



dt 



■ ,R V^(-70 + 28/3 2 +3/3 4 



(35) 



One obtains dv CR /dt > 0, contrary to what one would ex- 
pect. In Sec 4, cf. Eq (|40p . we show that that resonant diffu- 
sion has dominant effects on the evolution of CR streaming 
motion, which ensures dv CR /dt < 0. 



4 RESONANT DIFFUSION 

Waves generated from the nonresonant CR-streaming insta- 
bility can interact with CRs in resonance at p = x. Although 
the fraction of CRs in the resonance is small, we show that 
diffusion through such resonant interactions has a dominant 
effect on the CR streaming motion. Provided that the phase 
speed of the wave is <C cfi, resonant diffusion is equivalent 
to pitch angle scattering, in which particles change their di- 
rection of motion without gaining or losing energy. Since the 
nonresonant instability is in the regime i<C 1, resonant scat- 
tering involves mainly CRs with pitch angles a ~ 7r/2 ± x. 
Since one deals with purely growing waves that can be re- 
garded as a quasi-mode with ui — > 0, the assumption of low 
phase speed <C c/i is valid even for particles with pitch an- 
gles near 7r/2. The standard formalism involves writin g down 
the d iffusion equation involving pitch angles only (I Melrose! 
1986). Using F(u,n) = F(u ]} , u± ) with the replacements 
U|| = ufj, and 111 = tt-y/ 1 — /i 2 and retaining only the pitch 
angle scattering part in (|A12|) . one obtains 



dF(u, jj) 



dt 



d_ 

dfx 



F{u,n), 



D = -- 



7T n SB 2 
4 7 "b 5 "' 



where 

5Bl _ 
Mo 



dkiki 8Bi 



(2tt) 2 ho 



(36) 



(37) 



is the density of magnetic energy at the resonant wavenum- 
ber k r = Q/cu\/j,\. Using ([1]), it can be shown that the res- 
onant diffusion does not contribute to change in (u 2 ) and 
(uj_). Thus, the rate of change in parallel and perpendicular 
pressures is determined by nonresonant diffusion. 
The rate of change in the streaming speed is 

% = -Dv CR , (38) 
where 



D 



dfi(l 



H 2 )D: 



71 

T 



7 B 2 k r g0 ' 



(39) 



where 5B r is set the value at /1 ~ 7r/2. The approximation 
in is obtained by writing SB 2 k = (2Tv) 3 VSB 2 8(k ± )5(k r - 
ko)/k±, i.e. one assumes that the magnetic fluctuations are 
centered at the wave number ko- Here r 9 o is the gyroradius 
of CRs with an average Lorentz factor 70 = (7). 

Since the rate of the change in streaming motion due to 
resonant diffusion is (dv CR /dt) r ocil/y and the correspond- 
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ing rate due to nonresonant diffusion is (dv CR /dt) nr oc r/7 2 
one has 



(dv cn /dt) r 



(dv cn /dt) nr 



5 x 10" 



3 7« 



SB 2 
JW 2 



1 



korgo 



(40) 



Here SB is the typical amplitude of the magnetic fluctua- 
tions generated by the nonresonant instability. The right- 
hand side is larger than 1 for kor g o > 20(F/ 7 of2), where one 
assumes 5B r ~ SB and B,/B ~ fcor 9 o- This condition is 
satisfied for the typical parameters relevant for SN shocks 
(cf. Sec 5). If one writes SB oc exp(Fi) and takes the limit 
t > Eq J38j can be solved to yield 



SB 



12 r 7 /« cr0 

"5F ~ — T7 fc or 9 oln -2S_ 
iJ 7r \ u-,,. 



(41) 



5 MAGNETIC FIELDS IN SNRS 

There are a number of SNRs whose post-shock magnetic 
fields have been estimated from X-ray observations and 
found to be in the range ^ 10 -8 T, m uch higher than 
the typical value in the ISM (~ 10~ 10 T) (IVink fc Lamind 
120031 ; IVolk, Berezhko fc Ksenofontovl 12009 ) . Although one 
cannot rule out the possibility that such strong field 
is due to compression of strong pre-existing magnetic 
fields in the stellar wind of the progenitor star, the 
observations favor the interpretation that the magnetic 
fields are due to an instability at the shock. Obser- 
vations show that the magnetic field strength s t rongly 



vations snow mat tne magnetic neia strengtn s trongly 
depends on the shock speed (Vink fc Lamind 120031 ; 
IVolk. Berezhko fc Ksenofontovll2005 ). The post-shock mag- 
netic energy density in SNRs whose shock speeds have 
been estimated show a tendency SB 2 oc vj, where v B is 
the s hock speed (|Volk. Berezhko fc Ksenofontovll2005l ; IVinkl 
2008). The argument that the magnetic pressure is propor- 
tional to the ram pressu re (oc v 2 ) at the shock front suggests 
6 = 2, while iBelll (|2004l ) , s model predicts 8 = 3. So far, cur- 
rent observations cannot differentiate these two because of 
the nar row dynam ic range of the observed properties of SN 
shocks (|Vinkj|2008l V 



5.1 SN shocks 

We consider a forward shock traveling at a velocity v s . The 
shock undergoes initially a free expansion phase and after 
sweeping up sufficient mass, it enters a Sedov phase. In the 
Sedov phase, the velocity is given by v s — 2R s /5t, where R s 
is the shock radius. The total kinetic energy of the shock is 
E = (4/3)Rspov 2 , where p is the density of the ISM. On 
eliminating R s in favour of t, the shock velocity becomes 



1/5 



\125p J 



(42) 



One can show that development of the nonresonant instabil- 
ity is faster in the early Sedov phase. The maximum growth 
rate is oc J C r , which can be written in terms of the ram pres- 
sure in the up-stream region. Thus, one expects the growth 
rate to be a strong function of the shock speed. Since the 
ram pressure is ~ pov 2 , the CR pressure can be expressed 
as Pcr = Vapov 2 , where r] s ^ 1 is the acceleration efficiency. 



We consider precursor CRs with the initial streaming veloc- 
ity v CR0 ~ v s - Since Pcr = ncRmc{uv) /3, one has Jcr = 
3rj 3 qp v 2 vcR/(mc{uv)) w 3r} a qpoVsVCR/[mcui ln(M 2 /ui)], 
where (uv) « cuq and uq is the average momentum given by 



f , , p-1 1 - {u 2 / Ui ) 2 -p 
ug{u)du = ui - — — — . 

p — 2 1 - (U2/M1) 1 p 



(43) 



The average Lorentz factor of CRs is 70 = (1 + u 2 ) 1 ' 2 . For 
p = 2 and 112 > m, one has 110 = ui \n(u2/ui). Therefore, 
the faster the shock speed the higher the growth rate. 



5.2 Saturated magnetic field vs shock speed 

One may estimate the maximum magnetic field, SB m , of the 
wave by equating the magnetic pressure to the CR pressure, 
SB^/2p = Pcr- This leads to SB m w (2poVsPo) 1/2 v s , or 



8B m w 3 x 10~V /2 



Vl0 6 m- 3 / V 



5 x 10 6 ms- 1 



T.(44) 



Eq (|44|) is a rather optimistic estimate and calculation of 
the saturation due to the feedback effects of the instabil- 
ity on the streaming anisotropy generally leads to a lower 
saturation level than (|44[1 . 

One can calculate the saturated magnetic field due to 
the diminishing of the streaming anisotropy from (|41|l . This 
involves calculation of the cut-off streaming speed, denoted 
by v' CR , at which the instability turns off. From (|18[) . one 
obtains 



5c U B 
24rj s p v1 



k\\r g0 , 



(45) 



where Ub = B 2 /2po is the energy density of the mean mag- 
netic field. At a given parallel wavenumber feii , wave growth 
stops when the streaming speed is reduced to v CR ^ v' cn . 
One may assume v CR0 ~ v s . Eq (|45[) gives 



5 x 10 l 



k\\r g0 



. 5 x 10 6 m t 



itr 



'T 



B 



( no ) 
VlO 6 m- 3 / 



.10-m-"' (46) 
For convenience, one assumes here that CRs and ions in the 
ISM are all protons. 

The growth rate (|17p is estimated to be V w 
(6r? s fc r 90 ) 1/2 (i; s /c) 3/2 n/7o. From gTJ one obtains 



SB 2 
B 2 



v^) 1/2 (5) 3/2 (W /2 



hi 



«cr0 



(47) 



Eq (|47p implies that significant amplification 5B 2 /B 2 S> 
1 is possible if one assumes a large fcor s o- However, one 
should point out that the quasilinear theory is based on the 
weak turbulence assumption that requires the relevant wave 
fluctuctuations be small. Nonetheless, since the saturation 
mechanism discussed is quite generic, i.e. the saturation is 
attributed to the reduction in CR streaming, one expects 
that extrapolation of Eq (|47)> to the 8B 2 /B 2 > 1 regime 
would still provide a reasonable estimate for the saturation. 
To estimate the upper limit to the saturation one assumes 
that the Larmor radius is limited by the scale length of the 
acceleration region, which is assumed to be the shock's ra- 
dius, R s - The characteristic wavenumber ko is limited by the 
condition (JTSj) , which gives ko ~ (w p ,cr/c)(« cr /ua), where 
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u Pi cr is the plasma frequency of CRs. This corresponds to 



(fcor 9 o)max 

and fc r 9 o 



10" - W. For r) s 



0.1, v s 



5 x 10 6 ms 



10", one has SB 2 /B 2 ~ 39. The saturated 



magnetic field can exceed the background field by a modest 
factor. If the CR pressure is a fixed fraction of the ram pres- 



sure, one has n r 



. 7o oc v s 



Since the shock radius can be ex- 



pressed as R s = (3-E/4 / oof?) 1/ ' 3 , one has koR s oc vt^ , where 



1/3 



feo oc v c 



■ Eq (|47[) implies that the energy density of the 



saturated magnetic field increases with the shock speed as 
SB 2 oc v 2 . Remarkably, this predicts the same shock-speed 
scaling as the one from the equipartition argument, i.e. the 
magnetic pressure equal to the CR pressure (cf. Eg |4"4")| . 



6 CONCLUSIONS AND DISCUSSION 

We dev elop a kine tic version of the CR streaming insta- 
bility of iBelll i|2004h , who treated it using MHD. We show 
that the nonresonant growth is not due to the CRs them- 
selves, but rather to a current in the background plasma 
that neutralizes the current due to streaming CRs. Satu- 
ration of the CR streaming instability is discussed in the 
quasilinear theory. Growth of the nonresonant instability 
leads to particle diffusion that suppresses the instability. 
In the quasilinear diffusion formalism, the saturated mag- 
netic field can be calculated analytically. We consider both 
resonant and nonresonant diffusion, with the latter being 
treate d in a similar way to that us e d for the fireho se insta- 
bility l|Shapiro fe Shevchenko|[l9ol lDavidson|[l972T ). In our 
model the CR pressure is set to a fixed fraction of the ram 
pressure at the shock front. It is shown that evolution of the 
parallel and perpendicular pressures is due to nonresonant 
diffusion. The direction of the evolution is that the parallel 
energy increases and perpendicular energy decreases. This 
is in contrast to the firehose instability which causes the 
perpendicular pressure to increase and parallel pressure to 
decrease. One should note the difference in the sources of 
free energy that drive the instabilities; the free energy for 
the firehose instability is an excess of parallel pressure while 
the free energy for the nonresonant instability considered 
here is streaming motion of CRs. It is shown that satura- 
tion is determined by resonant diffusion which reduces the 
streaming motion. The saturated magnetic fluctuations can 
exceed the background field provided that the shocks are 
sufficiently fast and saturation occurs at short wavelength 
1/fto < I- 

There have been numerical simulations of the streaming 
instability in both the MHD approximation and kinetic 
theory, which confirm rapid growth of the instability under 
the physical conditions of young SN shocks. However, these 
simulations predict different saturation levels. One of the 
main difficulties in numerical simulation of the instability 
is the huge difference between the CR number density 
and background plasma density, with n CR /no ~ 10~ 5 . 
A particular approximation is usually adopted to make 
numerical simulation of such system feasible. For e x- 
ample, in iNiemiec. Pohl. Stroman. fe Nishikawal j2008)'s 
numerical model, the density ratio is set artificially to 
a large value. Because the existing numerical models 
use different approximations, and it is unrealistic to 
make meaninful comparison among their results. MHD 
simulations generally predict a saturation at strong mag- 



netic fields (|Belll |2004| ; IZirakashvili. Ptuskin fe VoUd 
2008), while kinetic particle-in-cell (PIC) simula- 
tions appear to predict much lower saturated mag- 
netic fields ([Niemiec. Pohl. Stroman. fc Nishikawal 120081 ; 
iRiquelme fc Spitkovskvl " 20081 ). The analytical calculation 
presented here predicts a relatively modest amplification 
that gives rise to magnetic field fluctuations similar to 
the level predicted by the recent kinetic simulations 
(|Riquelme fc Spitkovsk"vll2008h . 
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APPENDIX A: QUASILINEAR DIFFUSION 

Using the similar method discussed in iDavidsonl l|l972t) . we 
write the distribution function as / ~ /'°' + p 1 ' with 
a linear function of electric and magnetic fluctuations, 8E 
and SB, of the wave. A formal approach to determine the 
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feedback effect of the wave on the particle distribution in- 
volves calculation of time-dependence of the distribution av- 
eraged over the random phase, denoted by F = (f )r, 
where {■ ■ -}r is the random phase average. 



Al General formalism 

From the Vlasov equation one obtains 



dF /, crn c „, df (1) 



0, 



(Al) 



where d/dt — d/dt + vo ■ d/dx, «u is the drift velocity. Let 
£ = (SE,5B,f^) and the corresponding spatial Fourier 
transform be fjj. = (SE k ,5B k , fj, ), defined by 

£(*,*) = / j^£k{t)exp(ik-x^, (A2) 
where we assume WKB for the time dependence 



Cfc (t) oc exp 



-iipk(i) 



u) h (t)+iT(t) 



dt'. (A3) 



Here uj k is the wave frequency and Y is the growth rate. One 
adopts the convention u- k = — u k , £,- k = 

In cylindrical coordinates, the velocity can be expressed 
as v = (v± cos <f), v± sin <j>, v«), where <f> is the gyrophase. The 
3-axis is assumed to be along the magnetic field. For a trans- 
verse wave propagating parallel to the magnetic field, one 
has SE± = SEkx ± SE ky /i = ^f(iu)/k\\)(dB k x ± SB ky /i) = 
^(iuj/k\\)8B±. In the limit k±v±/£l <C 1, the perturbation 
of the distribution can be written as 



I f+e^ + f-e 



f (i) 

Jk ~ 2 



From the linearized Vlasov equation one obtains 
±qSB ± 



(A4) 



f± 



UJ + iF — kuv« ± fi/7 




d P\\ 



lo + iV\ d 
" fc || / dp± 



F. 



Substituting (|A4I) and (|A5|) into (| Al|) . one derives 



dF q I dk 

~dt = ~4iV J (2tt) 3 



_9_ 

a P|| 



w — ir \ i 9 



— ^ — P-L 



(5B* + f + -5Blf- 



where one uses 

<<5B fc ■ 8B k ,) R = M!j5B fe | 2 5(fe + ft') 

With \5B k \ 2 = |<SB 4 
dF q f dk 



(A6) 



(A7) 



|<5-B_j 2 , one obtains 



tft 2iV y (2tt) 3 



x 



w-iT\ 1 9 
fcii I p± dp± 



u> + iF — k\\Vn 



(w + ir-fc||U||) 2 -fi 2 /7 2 
d 



d P\\ 



to + iV\ d 
ku J dp± 



(A8) 



For practical purposes, we set ui ~ here. 
A2 Nonresonant diffusion 

The quasilinear diffusion equation in the long wavelength 
approximation can be obtained using the following expan- 
sion in I/O.: 

u + iF — k\\V\\ k\\V\\ iF 

DP' 



(u + iF- k }l v {] ) 2 -Q 2 /-y 2 



Q 2 



(A9) 



with Q = Q/y. We consider a more general expansion in 
r/|fe|,«n - f2| > 1: 

iF — k«v« 



1 



(iF- fc||?;||) 2 - ft 2 



A) II || f^*^ 
x I iF — kiivu 



2irfc 2 « 2 

k 2 V 2 - Ci 2 



(A10) 



Integration on the right-hand side of (|A8|) includes terms 
oc fey or oc 1/fey . For turbulent spectra with symmetry ku — > 
—fen, one has 

dfcy (terms oc fey or oc l/fe||) = 0, (All) 
and only terms oc iF remain. 

A3 Resonant diffusion 

The standard procedure to write down the resonant diffusion 
involves setting F = 0. Eq (|A8fl can be expressed in the 
following familiar form 



^ ~dt = ~ 



dF d I n d _ d 



d_ 

d P\\ 







1 

p± dp 



D i 



dpi 



<9p_l 



(A12) 



where are the diffusion coefficients. 



